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Temperature correction to the Casimir force is considered for real metals at low temperatures. 
With the temperature decrease the mean free path for electrons becomes larger than the field 
penetration depth. In this condition description of metals with the impedance of anomalous skin 
effect is shown to be more appropriate than with the permittivity. The effect is crucial for the 
temperature correction. It is demonstrated that in the zero frequency limit the reflection coefficients 
should coincide with those of ideal metal if we demand the entropy to be zero at T = 0. All the 
other prescriptions discussed in the literature for the n = term in the Lifshitz formula give negative 
entropy. It is shown that the temperature correction in the region of anomalous skin effect is not 
suppressed as it happens in the plasma model. This correction will be important in the future 
■ cryogenic measurements of the Casimir force. 
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> . I. INTRODUCTION 

in 

■ Attraction between parallel metallic plates predicted by Casimir in 1948 Q| (see for a recent review) has been 
measured in recent experiments |^5§i^7§|j7|J^ with high precision. To calculate the force with the same precision, one has 

■ to take into account different corrections to the original Casimir force The corrections due to finite conductivity 
' of the plates and roughness of their surfaces can be large but there is no principal problems with them and the 

discussion is going around the material properties [p|l0|,ll,12 
f~i ' the finite temperature raised debate in the literature 1 13. 17, 18 
Qh, different authors has not been reached yet. 

^li • In the early period the temperature correction was found for ideal metal. It has been discovered that the correction 
^ following from the Lifshitz theory p6 1 of electromagnetic fluctuations in nonhomogeneous media did not agree with 
that found with different methods |27j2q ]. The problem originates from the n = term in the Lifshitz formula giving 
the classical contribution in the Casimir force from long wavelength fluctuations. Special prescription for this term 
was proposed by Schwinger, DeRaad, and Milton to reconcile different approaches. Namely, one has to take the 
limit of infinite permittivity before allowing the frequency go to zero. 

The problem arose again when researches tried to find this correction for real metals. It has been realized that the 
5J] plasma and the Drude models for the metal dielectric function gave different corrections, which did not coincide when 
\ the Drude relaxation frequency lOt is going to zero. Direct application of the Lifshitz formula with the Drude dielectric 
function gives large linear in temperature correction p6[ |, which contradicts to the torsion pendulum experiment 
Negligible for the experiments correction was found using the plasma model ]T^JT^ . However, the Drude model works 
much better for real metals and there is no a satisfactory reason to justify the use of the plasma model. 

It became obvious that something wrong with the n — Q term in the Lifshitz formula and while a solid theoretical 
reason for modification of this term is not found one has to use a reasonable prescription to define the n = term. It 
was proposed ||l^ to apply for real metals the same prescription as for ideal metal ||2^]. It was motivated by the fact 
that in the static limit the boundary conditions on the metal surface did not depend on any material property. Later 
it has been shown [ p^ that the classical term can be constructed from very general principles using the dimensional 
analysis, low frequency limit for the permittivity, and the known result for ideal metal. The same conclusion was 
made in a recent paper p5| ], where the authors stressed difference between boundary conditions on dielectrics and 
metals in the zero frequency limit. In this case the linear in temperature correction survives but it is suppressed by 
small additional factor. The correction is unobservable p3| in the conditions of torsion pendulum experiment |H] but 
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3 times larger than the experimental errors in the atomic force microscope (AFM) experiment Note that it 
makes agreement between theory and experiment better. 

Quite different prescription was proposed for real metals by Klimchitskaya and Mostepanenko ||2^. The driving 
idea was to modify the reflection coefficients in the n — term in such a way to make them coincide with that of 
plasma model in the limit lUt —^ 0. The reason for that was based on the property of the reflection coefficient for 
perpendicular polarization ri((^,q) considered as a function of two variables: continuous imaginary frequency C and 
the absolute value of the momentum along the plate q. It has integrated discontinuity in the point C = 0, q = 
which authors [po[ consider as unacceptable. The resulting temperature dependence does not include the linear in 
temperature term and negligible in conditions of the torsion pendulum and AFM experiments. 

Thus, at the moment there are three different approaches to the problem of temperature dependence of the Casimir 
force. Recently an interesting proposition to check the Nernst heat theorem for different prescriptions discussed 
in the literature was made. According to the third law of thermodynamics the entropy S of the closed system in 
equilibrium must go to zero in the limit T —^ 0. Based on the analytical result for the plasma model and numerical 
calculations for the Drude model a conclusion was made that the prescription of Schwinger, Deraad, and Milton 
cannot be applied to real metals because the entropy still finite in the T — > limit. On the contrary, it was concluded 
that for the prescription proposed by Klimchitskaya and Mostepanenko S ^ when T ^ 0. 

In this paper we intended to demonstrate that the real situation in the low temperature limit is different. The reason 
for this is that with the temperature decrease the mean free path for electrons I increases and at some temperature it 
inevitably becomes larger than the field penetration depth S. It is the range of anomalous skin effect when the space 
dispersion becomes important. Taking this effect into account modifies the temperature correction, which does not 
coincide with that found in Ref . ||30|] . Independent interest to the correction at low temperatures is stimulated by the 
progress in experiments which are going to cryogenic temperatures to improve the precision ||^ . 

The paper is organized as follows. In Section || we give definition of the Casimir free energy and explain the 
problem with the n = term. This term is written out for all prescriptions discussed in the literature. In Section |l| 
it is demonstrated that neither plasma nor Drude models can be used for real metals at low temperatures. The way 
to calculate the temperature correction using the impedance of anomalous skin effect is described. In Section ^ the 
actual calculations are presented. The results are applied to analyse the entropy behavior at T ^ and to find the 
value of temperature correction in cryogenic range in Section Our conclusions are collected in the last section. 



II. CASIMIR FREE ENERGY 



We start from the Lifshitz expression for the force per unit area between two parallel plates separated by the 
distance a at temperature T ||3l|. Simple transformation of this formula allows to write the free energy in the form 



n=0 / 



ri^ra], (1) 



where the prime over the summation sign indicates that the first n — term has to be taken with the factor 1/2 and 

^„ = nr, r = kT^fj = ^ = fu^a- (2) 

Here ^„ are the dimensionless Matsubara frequencies ^„ = Cn/^^a, where C„ = 2TTkTn/h. In Eq.(^ 2 iCn, o) are the 
reflection coefficients for orthogonal (ri) or parallel (r2) polarizations. These coefficients depend on the material via 
the dielectric function e (iCn) at imaginary frequencies. In the Drude model e (iCn) has the form 

£(»C)- ^ (3) 

C(C + UJr) 

and is defined by two parameters which are the plasma frequency ujp and the relaxation frequency u-r- This model 
nicely fits good metals in the whole frequency range excluding the interband absorption region to ^ ujp, which does 
not give considerable contribution especially at low temperatures. 

Consider first the problematic n — term. Expression was found by solving the equation for the Green function 



in nonhomogeneous media at finite temperature |31|. The explicit form of the refiection coefficients 2 depends on 



the boundary conditions for the Green function at the metal-vacuum interface. Note that these conditions one can 
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put separately for any Matsubara components. Typically the continuity of tangential components of electric and 
magnetic fields at the interface is demanded. The Drude model gives for the amplitudes ri_2 

where the function is defined as 



In the n = term we have to put = and the reflection coefficients are 

r2(0,y) = 0, r2(0,y) = l. (6) 

while for ideal metal both of these coefficients are equal to unit r\ = r\ = \. There is no way to reconcile the value 
of r\ with that of ideal metal since it does not depend on any parameters. This is why one needs to introduce a 
prescription for the n = term. The coefficients were used to calculate the temperature correction to the Casimir 
force in Ref. p^ . 

It was proposed to use the extension of the Schwinger, DeRaad, and Milton prescription to real metals. Since 
in the static field limit the boundary condition on the metal surface does not depend on the characteristics of a 
particular metal, the reflection coefficients were changed with that of ideal metal 

r2(0,2;)^l, r2(0,y) = l. (7) 

Different prescription was introduced by Klimchitskaya and Mostepanenko pO| ] who proposed to modify r\ in such a 
way that in the limit t^T- — > it coincided with the plasma model result 



r\(SS.y)^rH^y,y)^\^l^^2^\ , r^(0,g) = l. (8) 

V {y) + y^ + y) 

The motivation of (Q) and (||) was briefly described above and we are not going to the details referring to the original 
papers @|0|. 

Let us denote the n = term in (|l|) as J^o{a,T). Calculating the corresponding integral with the reflection 
coefficients (0-0) one gets 



kT 

J^oia,T) = -a^Ci3), (9) 



where Ci^) is the zeta function. The coefficient a has the following values for the prescriptions (||), (Q), and (H), 
respectively: 



(i) 1 (ii) 



a(-' = 1 - 4^ - V O f4V (10) 



ujp ujp({3) \uj, 
The function I2 {x) varies slowly with x and is given by Eq.(68) l2 



III. ANOMALOUS SKIN EFFECT 

Let us consider now the n > 1 terms in (0) when the temperature is going down. The plasma frequency ujp is defined 
by the electron density in metal and practically does not depend on temperature. On the contrary, the relaxation 
frequency u-r changes with T significantly. If T < To , where is the Debye temperature for a given metal, then the 
dependence uJr (T) can be presented as (see, for example, [p2[) 
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UJr (T) = UJr (0) + C,r2 + CphT^' 



(11) 



with the parameters ^^^(O), Cg, and Cph- Here the first term is connected with the scattering on lattice irregularities 
and impurities, the second one describes scattering on electrons, and the third term corresponds to the scattering 
on phonons. The first and second terms can dominate only at very low temperatures. The first term needs special 
discussion. The relaxation frequency is proportional to the material resistivity which disappears at T — > for perfect 
monocrystals. If we are going to check the third law of thermodynamics for our system, we have to take this equilibrium 
state and choose Wr(0) = 0. To describe the temperature behavior of real material used in the experiment, it can be 
wrong. For evaporated metallic films the residual resistivity can be significant due to large density of defects but it 
is much smaller than the resistivity at room temperature. 

One can easily see the result of Ref. without any calculations. At low temperature uj^l^a in (^) is going to zero 
faster than ^„ ~ T. Therefore, for all n terms one can use the plasma model. The reflection coefficients in the 
n — Q term (^) were prescribed to reproduce the plasma model at Wr — > 0. Therefore, the plasma model completely 
describes the situation at low temperatures. It is well known [0| that the leading temperature correction in this 
model behaves as and the entropy will go to zero as S' ~ T^. Ai^ other prescription will not agree with the Nernst 
theorem since the n = term does not coincide with that for the plasma model. For the entropy becomes negative 
and for (|^) it is positive but both of them are finite at T = 0. No doubt that at low temperature and for equilibrium 
state the relaxation frequency becomes negligible. However, it is wrong to think that the reflection coefficients (0) 
for n ^ will be the same at low temperatures. This is because the mean free path for electrons increases with 
the temperature decrease according to the relation I — vp IoJt{T)^ where is the Fermi velocity. Atj-some small 
temperature inevitably the relation I ^ 5 will be fulfilled, where 5 — c/ujp is the field penetration depthu. Then the 
condition I ^ 5 is equivalent to the following 

ujr (T) < —ujp = n. (12) 
c 

The frequency fl = {vp/cjiVp is often used as a characteristic frequency of anomalous skin effect. 

Thus, the local connection between the current and electric field, which is true for I < d, is broken at low temperature 
and space dispersion becomes important. This is the range of anomalous skin effect when the dielectric function 
s{uj) cannot be used any more for description of the metal. Instead the interaction with the field is defined by the 
surface impedance Z{ui) |^ , ^ . The relation Z{uj) ~ e~^^'^{uj), which holds for the normal skin effect, is broken. 

When the metal is described by the impedance the boundary condition on its surface will be 

Et = ZH(HtXn), (13) 

where Et, Hj are the tangential components of electric and magnetic fields, n is the unit vector normal to the surface 
and directed inside of the metal. It holds true while the impedance is small. If this boundary condition is used for 
the Matsubara components (n > 1) of the Green function, then the reflection coefficients can be written as follows 

Impedance in the range of anomalous skin effect one can find solving the kinetic equation for not in equilibrium 
distribution function of electrons in the electromagnetic wave or using simple qualitative analysis [^^ . 
In the frequency range uj < H, (strong anomalous skin effect) the impedance is |33p3] 



Z{u.)=e-/^{^y\ (15) 



Here v = (3vf, where /3 is a factor ^ 1, which is defined by the structure of the Fermi surface. Note that the relaxation 
frequency falls out from the impedance at all. This is because the field interacts mainly with electrons, which are 
moving in the surface layer of thickness S. For this reason significant is the effective conductivity fe// = a6/l, which 



'^Strictly speaking 5 depends on frequency. We assumed that iv — ^ ujt (T) [n ^ 0). If it is not the case, we have to take 
S — {c/ujp)y^2ujT/uj but the qualitative conclusion will not change. 
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does not depend on ujr- Actually we are interesting in the impedance at imaginary frequencies. For the strong 
anomalous skin effect it will be 

This expression is true for ^„ < fl/tOa while the range of important n in the sum ([^) is given by the condition ^„ ~ 1. 
For small fl/cua the relation (|l^) cannot be used for all important n. However, as we will see below the temperature 
dependent part of free energy always can be described by (p^. 

The anomalous skin effect will change not only the temperature correction but the main temperature independent 
part of the Casimir force or free energy. To find this change one has to know the impedance in the whole range of 
It is very important for comparison with the experimental data in cryogenic range but out of the scope of this paper 
and will be discussed elsewhere. 



IV. CALCULATION OF THE FREE ENERGY 



Now we are ready to find the sum of all > 1 terms in (|T]). This sum can be transformed using the Abcl-Plana 
formula p7|| . To avoid inconvenient integrals we first change the summation index to m = n — 1. In this way the 
problematic coefRcient ri(0, y) will not appear in the formula as it usually happens The resulting expression is 
longer than usual but simpler for calculations. Separating the temperature independent part one finds 



T{a, T) - To{a, T) = T{a, 0) 



kT 
87ra2 



r(l) 



r(2; 



(2; 



Here the temperature independent term !F{a^ 0) is 

he 



00 00 



Jd^Jdy [(f + y) In (1 - rl y + e-^-«) + n ^ ra] 



(17) 



(18) 



and the integrals ijf^ (j — 1,2; k — 1, 2, 3) are defined as 



= J dyy\n[l-r]iT,y)e-y] , 



-'2 



i oo 

j dt j dyy In [l - (rt, y) i 



= 2/m J -5^^ J dy{y + irt) In [l - r] (r + zri, y + irt) e-y-"'\ . (19) 

r 

In the temperature independent term ( |l8|) the important range of continuous variable ^ is ^ ^ 1 and, as was mentioned 
above, the impedance ( |l6| ) cannot be used to cover all the range of ^. But this term is not interesting for us. In the 
temperature dependent terms ^ = rt ~ r as one can see from ([l9|). We can use (|l^) in the integrals ( [l9| ) if r < f2/wa 
or equivalently 27rfcr < TiKl. This relation together with T < Tu will be supposed to be true. 

In all the integrals ( |l9| ) the lower limit in y can be changed into zero. Really, the error connected with such a 
replacement will be of the order of Inr. It will be considered as small for analytical calculations but will be taken 
into account in the numerical procedure. Then the integrals /^^^ depending on the reflection coefficient r\ ( p^ ) will 
be defined only by the parameter 

(20) 
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Even for small t the typical value of this parameter is large A ^ 1. Only for very low temperature kT <C 
{u!a/uJp)^hfl/2Tr it becomes small {A 1). It is in contrast with the parameter B which defines the integrals 



r(2) 



V LO 



2 



1/3 



B-[-^r^] . (21) 



This parameter is always small in the range of interest t ^ 1. 



A. Free energy in the limit A 



Let us analyse now the analytical dependence of the integrals ( |l9[) from the parameters A, B at very low tempera- 
tures when A <C 1. We start from Neglecting the correction Inr this integral can be written 



oc oc 

I^i^ = j dyy\n{l - e"^) + j dyy\n 



(e^ - 1) 



(y + A) 



(22) 



The first integral here is equal — C(3). It describes contribution of the ideal metal. For the second integral one can 
find the leading terms in the limit A — s- 0. The important contribution gives the range y < 1 but the main one is 
collected near y ~ \/A. Then 7^^'' can be presented as 



1 

/««-C(3) + J dyyln(^l + ^y 



The result is the following 

/f) = -( (3) - 2A{ItiA + 2ln2 - 1) + O {A^) . (23) 

The integral /j^' can be estimated in the same way but one has to change A ^ A t^^^ and then integrate it over t. It 
gives 

4'^ = -C (3) ~lA(lnA + 2\n2-^^+0 {A') . (24) 

The same procedure can be applied for /g^'' but realization is more complicated. Contribution of the ideal metal 
(zero impedance) is negligible in our approximation since the leading term is known p7| , p8| to be t^. Correction due 
to nonzero impedance depends on A. In the same approximation as for ( |2^ ) and (|^ one has 



= 2Im 



dt 



1 



dy {y + iTt) In ( 1 + ^ ^ 



The inner integral is calculated analytically. The contribution from the term proportional to irt will be of the order 
of t\/A which is small in comparison with the leading terms. All the rest can be written as 



= 4A(gi \nA + ga) + O (rVl) 



(25) 



The constants qi^2 here are 



91 



dt(H-t2)i/^ sin I? 



92 = 



dt 1 + 



(26) 
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where tan Si? = t. Numerically these constants are qi = 0.0137, q2 = 0.0191. 

(2) 

The integrals /j; can be calculated quite similar. Since B is always much smaller than A one can completely 

(2) 

neglect the terms containing B and for Ij^ one has 



r(2) 



-C(3), /: 



(2) 



-C(3): 



7-(2) 



(27) 



Substituting (^-25, 2^ into ( |T7| ) and taking into account explicit expression ^ for To{a, T) one finds the temperature 
correction to the free energy in the limit A ^ Q 



AT (a, T)^T (a, T) - T (a, 0) 



kT 



(l-a)C(3)+^ 



In A 



In2-I 



4,2 



(28) 



This expression is true for very low temperature (see Eq.(^0|)) and is appropriate only to check the Nerst theorem. 
It will be discussed in the next section. For realistic low temperatures, which will be explored in the near future 
experiments, the typical value of A is large even for small r and we have to analyse the opposite limit A^ 1. 



B. Free energy in the limit 1 

To calculate the integrals /[.^^ in the temperature range when the parameter A is large but the temperature still 
small (r <^ 1), it will be convenient again to separate the contribution of ideal metal and neglect the correction 

~ Inr due to change of the low limit in integrals into zero. Then for one has the same representation ( p^ ) and 
similar for the other integrals. Because A is large one can expand the logarithm in series and collect the coefficient 
at different powers of 1/A. For I^^"^ this procedure can be done straightforward and we find for the first two terms 

/j^^ = -C(3) + 8C(3)(^i--^). (29) 

In case of /j^'' the role of A plays At^^^. In this integral after logarithm expansion one has to make first the 
integration over t not to run onto divergencies. The result is the following 



4'^ = -C(3) + 12C(3)(^^-^j. (30) 

In the third integral instead of A appears A{1 + ity^^. In our approximation one can completely neglect irt in 
comparison with A or with y. Making the same procedures one gets 

4^' = 16C(3)(^-5), (31) 

where the coefficients are defined as 

oc oo 

f , sin I? /■ , sin2i? 
Pi = / dt — , P2 = / dt — . 32 

J (e2'^*-l)(l+<2)l/6 J (g2^t_i) (1+^2)1/3 

Here the angle "i? is defined as in (p6|). Numerically the coefficients are pi = 0.0133, P2 = 0.0262. 

The integrals 7^, giving the contribution of parallel (r2) polarization depends on the parameter B which is small 
while T <C 1. Therefore, they give the same r esult (|7|) as in the case of small A. Collecting all together one finds the 
temperature correction to the free energy in the limit A 3> 1 



1 12 



kT 



(1 



1 - 2pi 15 - 12p2 



A A2 



(33) 



The temperature correction to the Casimir force can be easily found via AjF. Correction to the force between sphere 
and plate AFgp is directly proportional to AT if we adopt the proximity force theorem 

AFsp{a,T) = 2TTRAT{a,T), (34) 
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where R is the sphere radius. The temperature correction for the force between two plates A.Fpp is 



8 kT 
AF,,(a,r) = --A^(a,T)^^C(3) 



(l-«)-8(-^--2— ^ 



(35) 



Thus, we found the asymptotics for AJ^{a,T) in the limits A ^ ( p8|) and A ^ oo (|33|). Now let us calculate 
numerically this function in the transition region. 



C. Numerical calculation of the free energy 



For numerical calculation we parametrize the temperature correction to the free energy in the following way 



kT 

A.F (a, T) = ^ [(1 ~ a) C (3) - G {A, r)] 



(36) 



The absolute scale of the correction is given by the factor kT/Sira^ then G{A, t) can be treated as the relative 
temperature correction if a = 1. The function G{A,t) is calculated numerically as a function of A for a few given 
values of r. There is no need to separate B as an independent parameter since it is connected with the other two 
as -B = t'^/A. Alternatively at fixed material parameters ujp and v it can be found as a function of T at a fixed a. 
Actual calculation has been done using the relation 



GiA,r) = -l (/« + /f ) - 2C (3)) + (4^' + /f ) - 2C (3)) - (4^^ + . 

The integrals were calculated according to ( p^ ) with the absolute precision of 10^^. F ig.[l| shows the function G {A, 0) 
and its asymptotics a,t A —> and A ^ oo, which can be extracted from ( |2^ ) and (p3^. One can see that there is 
good agreement of numerical calculation and analytical asymptotics. If we consider nonzero but small r, the result 
will change only shghtly due to correction ~ Inr. For iVp = 1.37 • 10^^ rad/s and v = 1.5 ■ 10® cm/s (gold) G{a, T) 
as a function of T is shown in Fig.|| at a few values of a. This figure shows that the temperature correction is always 
significant in the interesting range of temperatures and distances between plates. 



V. DISCUSSION OF THE RESULTS 



A. Entropy in the limit T ^ 

We already saw that the anomalous skin effect changes behavior of the free energy at low temperatures. Now we 
are able to answer the question which of the prescriptions (||), (0) or (||) agrees with the third law of thermodynamics. 
In the low temperature limit entropy can be calculated from (EM) 



dT Sua 



(a - 1) C (3) -\a{C^+ Aq^ hi A + ln2 - ^ + Aq^ + <Zi 



(37) 



Parameter a corresponding to different prescriptions used for the n = term in the Lifshitz formula is given by Eq. (^^ . 
If this term is calculated as it appears in the Lifshitz formula without any modification, then a = a^*'* = 1/2. In 
this case the entropy is finite at T = (A = 0) and, moreover, it is negative. This conclusion coincides with that 
made in Ref. ||30tl , where at low temperatures the authors described the plate material with the plasma model. It can 
be considered as an additional physical argument that unmodified Lifshitz formula cannot be used to get temperature 
behavior of the Casimir force. The previous argument jl^,^ was that in this approach the n = term in (Q) does 
not depend on any material parameter (see (^ with a = aP^ ) and for this reason it cannot be reconciled with the 
ideal metal result, which is two times larger. These arguments show that the Lifshitz formula is really in trouble. It 
obviously has to be modified but after a few years of active discussion still there is no solid theoretical understanding 
what is wrong and how one can do this. 

Prescription proposed by Klimchitskaya and Mostepanenko pO| ] was designed to reproduce the plasma model result 
for the n = term in the limit lOt Q. This limit was supposed to be realized at low temperatures but, as was 
explained above, the material has to be described in this case by the impedance of anomalous skin effect rather than 
the plasma model permittivity. For this reason the entropy for a — a^"'^ is not going to zero when T 0. It is easy 
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to see from ( ^ ) and (|T^) that S is finite and negative at T = 0. Therefore, this prescription is also does not obey 
the Nernst theorem. At higher temperatures, when the Drude model is valid, the n = term given by (|9|,[l0|) depends 
separately from LUp and lUt- Physically it is unacceptable since at low frequencies the only parameter characterizing 
a metal in respect to electromagnetic field is the conductivity cr ^ i^fj^r | 3^ . 

If we use the prescription then a — a'") = 1. In this case the temperature independent term in is 

canceled completely and entropy is going to zero with the temperature as T^/'^lnT. It is the only prescription which 
is in agreement with the third law of thermodynamics. This prescription was criticized [[L7 U 20 I on the basis that the 
n = term does not depend on any material parameter (see (||)). We already explained [ll9| that in the static (long 
wavelength) limit the boundary condition for metals do not include a particular metal parameters. More specifically, 
from dimensional analysis it follows that a in (|p) can be a function of the only variable |l9[| 

a = a — TT- ■ (38) 
V / 

The characteristic frequency which appears here is huge Up/uJr ^ 10^^ rad/s (for good metals at room temperature). 
Of course, one can expand a in a series 



a « a (0) + a' (0) 

but the correction to a (0) becomes important only at microscopic distances between the plates a 10~® cm where 
the macroscopic Casimir force is not defined. The value of a(0) should coincide with that for ideal metal which 
is known to be 1. In this way the prescription (^ is reproduced without direct reference to the Drude reflection 
coefficients (^). 




B. Temperature correction to the Casimir force 



The results of previous section for the temperature correction to the free energy or equivalently to the force (see 
( p^ ) and (p5|)) show that for a = 1 the correction always increases the absolute value of the force (the correction is 
negative so as the attractive force). The correction is relatively large ~ kT/STra^. It is unlike to the plasma model 
where additional small factors or Wa/wp fl^ appear if one uses prescriptions (||) or respectively. In this 

sense one can say that in the range of anomalous skin effect the temperature correction becomes large. The relative 
correction G(a,T) in ( ^6|) is maximal (w 0.53) at the temperature 

fcT™«18^f^4V (39) 



In the ideal metal limit — > oo this temperature is going to zero and instead of large correction appears usual ideal 
metal correction |^,^ given by 

G(a,T).((3)(^)'-^!(^)'^ (40) 

Let us describe the temperature range where the result of this work will be applicable. Consider first the condition 
I 3> S which guarantees that anomalous skin effect plays important role. It can give us the upper limit on the 
temperature. We approximate lOt{T) with the Bloch-Griineisen formula 



u:r{T) _ F^{T/Td) 



iJriTo) \ToJ F5(To/Td) 



(41) 



where Tq is some fixed temperature, for example, Tq = 0° C. This formula does not take into account scattering on 
the defects and electrons, which can be important at very low temperatures, but it is good to find the upper limit 
on T. The value of ujriTo) can be fixed via the material resistivity p using the relation ujr = eoui^p, where Sq is 
the permittivity of vacuum. For gold parameters LOp = 1.37 • 10^^ rad/s, vp = 1.4 • 10^ cm/ s, p{To) = 2.06 pil ■ cm 
we found l/S > 5 for T < 113° K and l/S > 10 for T < 67° K. The other condition 2nkT < Un ensures that a 
specific expression for the impedance ( p^ ) will be true. It restricts the temperature by the value T < 77.5° K (gold). 
Therefore, already for liquid nitrogen our result for the temperature correction is applicable. 
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VI. CONCLUSION 



We considered the temperature correction to the Casimir free energy (force) in the low temperature range. The 
aims of this analysis were to investigate the behavior of entropy in the limit T ^ and find out if it is possible to 
neglect the temperature correction in the future low temperature experiments. The main observation of this work 
is that at low temperature T < Td (Debye temperature) neither plasma nor Drude models for material permittivity 
are good for real metals. The reason is that with the temperature decrease the mean free path for electrons becomes 
larger than the penetration depth for electromagnetic field in metal. This relation between parameters is realized for 
the anomalous skin effect for which description of metals with impedance is more appropriate than with the dielectric 
function. This change in the description is important for the temperature correction. 

It is known that the first n — term in the Lifshitz formula for the Casimir force is controversial. A lot of discussion 
in the literature is going on around this term. The problem is how one can define this term correctly. At the moment 
there are three different approaches corresponding to three different results for the correction. It was proposed to 
use the third law of thermodynamics to choose one of the approach ]30| ]. We analysed the entropy behavior in the 
limit T — !■ for discussed in the literature prescriptions and came to a conclusion that one has to define the reflection 
coefficients in the static field limit as for the ideal metal rf ^ 1, = 1 to get agreement with the third law of 
thermodynamics. It is in contrast with the conclusion of Ref. |3^], where the plasma model was used for metals at 
low temperatures. 

To reduce the noise, the future experiments on precise measurement of the Casimir force will explore the low 
temperature range. Previously metals at low temperatures were described by the plasma model because the Drude 
relaxation frequency decreases fast with temperature. The plasma model predicts a negligible temperature correction 
to the force. We demonstrated that the anomalous skin effect makes drastic change. In all interesting from the 
experimental point of view range of temperatures and distances between bodies the correction is not negligible on the 
precision level of modern experiments. 

We do not consider our result as final solution of the problem with the temperature correction. This is because a 
real solution should not be based on any prescription. The problem indicates that something wrong with the Lifshitz 
formula and efforts have to be directed on the careful analysis of this formula. It becomes more and more clear that 
the essence of the problem lies in the boundary conditions on the metal surface. The impedance condition ( [T^ ) applied 
to the n = term reproduced the right result for the reffection coefficients = = 1 |2^. However, it is not clear 
why there is no a smooth transition between continuity of tangential components of E and H and the impedance 
boundary condition. 
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